Relationship between endogenous hormone levels of grapevine callus cultures and their morphogenetic behaviour

Dark callus cultures from leaves and anthers of three grapevine genotypes (Seyval blanc, Riesling and Trollinger) were propagated in vitro for almost two years on medium deprived of growth regulators. Three different callus lines originated from these initial callus cultures: A dark non-embryogenic one, another composed by clusters of somatic embryos and a friable, whitish callus line that can produce somatic embryos under appropriate circumstances, Endogenous hormone levels viz, indole-3-acetic acid (IAA), abscisic acid (ABA), gibberellins A1, A3 and A20, zeatin/zeatin riboside and N6(∆2-isopentenyl) adenine/ N6(∆2-isopentenyl) adenosine, were analysed in the different callus types. Only higher ABA levels correlated with the morphogenetic capacity of the cultures. When endogenous hormone levels were analysed in the line constituted by conglomerates of somatic embryos at different dates during the whole culture period, higher levels of IAA and ABA were found during the winter months even if the cultures were maintained permanently under constant temperature and photoperiod conditions. A 4-week chilling period led to a reduction of the endogenous ABA level

Monolithic simulation of convection-coupled phase-change - verification and reproducibility

Phase interfaces in melting and solidification processes are strongly affected by the presence of convection in the liquid. One way of modeling their transient evolution is to couple an incompressible flow model to an energy balance in enthalpy formulation. Two strong nonlinearities arise, which account for the viscosity variation between phases and the latent heat of fusion at the phase interface. The resulting coupled system of PDE's can be solved by a single-domain semi-phase-field, variable viscosity, finite element method with monolithic system coupling and global Newton linearization. A robust computational model for realistic phase-change regimes furthermore requires a flexible implementation based on sophisticated mesh adaptivity. In this article, we present first steps towards implementing such a computational model into a simulation tool which we call Phaseflow. Phaseflow utilizes the finite element software FEniCS, which includes a dual-weighted residual method for goal-oriented adaptive mesh refinement. Phaseflow is an open-source, dimension-independent implementation that, upon an appropriate parameter choice, reduces to classical benchmark situations including the lid-driven cavity and the Stefan problem. We present and discuss numerical results for these, an octadecane PCM convection-coupled melting benchmark, and a preliminary 3D convection-coupled melting example, demonstrating the flexible implementation. Though being preliminary, the latter is, to our knowledge, the first published 3D result for this method. In our work, we especially emphasize reproducibility and provide an easy-to-use portable software container using Docker.Comment: 20 pages, 8 figure

FEniCS-HPC: Automated predictive high-performance finite element computing with applications in aerodynamics

Developing multiphysics finite element methods (FEM) and scalable HPC implementations can be very challenging in terms of software complexity and performance, even more so with the addition of goal-oriented adaptive mesh refinement. To manage the complexity we in this work present general adaptive stabilized methods with automated implementation in the FEniCS-HPC automated open source software framework. This allows taking the weak form of a partial differential equation (PDE) as input in near-mathematical notation and automatically generating the low-level implementation source code and auxiliary equations and quantities necessary for the adaptivity. We demonstrate new optimal strong scaling results for the whole adaptive framework applied to turbulent flow on massively parallel architectures down to 25000 vertices per core with ca. 5000 cores with the MPI-based PETSc backend and for assembly down to 500 vertices per core with ca. 20000 cores with the PGAS-based JANPACK backend. As a demonstration of the power of the combination of the scalability together with the adaptive methodology allowing prediction of gross quantities in turbulent flow we present an application in aerodynamics of a full DLR-F11 aircraft in connection with the HiLift-PW2 benchmarking workshop with good match to experiments

Dimension reduction in heterogeneous parametric spaces with application to naval engineering shape design problems

We present the results of the first application in the naval architecture field of a methodology based on active subspaces properties for parameter space reduction. The physical problem considered is the one of the simulation of the hydrodynamic flow past the hull of a ship advancing in calm water. Such problem is extremely relevant at the preliminary stages of the ship design, when several flow simulations are typically carried out by the engineers to assess the dependence of the hull total resistance on the geometrical parameters of the hull, and others related with flows and hull properties. Given the high number of geometric and physical parameters which might affect the total ship drag, the main idea of this work is to employ the active subspaces properties to identify possible lower dimensional structures in the parameter space. Thus, a fully automated procedure has been implemented to produce several small shape perturbations of an original hull CAD geometry, in order to exploit the resulting shapes and to run high fidelity flow simulations with different structural and physical parameters as well, and then collect data for the active subspaces analysis. The free form deformation procedure used to morph the hull shapes, the high fidelity solver based on potential flow theory with fully nonlinear free surface treatment, and the active subspaces analysis tool employed in this work have all been developed and integrated within SISSA mathLab as open source tools. The contribution will also discuss several details of the implementation of such tools, as well as the results of their application to the selected target engineering problem

Enhanced error estimator based on a nearly equilibrated moving least squares recovery technique for FEM and XFEM

In this paper a new technique aimed to obtain accurate estimates of the error in energy norm using a moving least squares (MLS) recovery-based procedure is presented. We explore the capabilities of a recovery technique based on an enhanced MLS fitting, which directly provides continuous interpolated fields, to obtain estimates of the error in energy norm as an alternative to the superconvergent patch recovery (SPR). Boundary equilibrium is enforced using a nearest point approach that modifies the MLS functional. Lagrange multipliers are used to impose a nearly exact satisfaction of the internal equilibrium equation. The numerical results show the high accuracy of the proposed error estimator

Investigation of wing crack formation with a combined phase-field and experimental approach

Fractures that propagate off of weak slip planes are known as wing cracks and often play important roles in both tectonic deformation and fluid flow across reservoir seals. Previous numerical models have produced the basic kinematics of wing crack openings but generally have not been able to capture fracture geometries seen in nature. Here we present both a phase-field modeling approach and a physical experiment using gelatin for a wing crack formation. By treating the fracture surfaces as diffusive zones instead of as discontinuities, the phase-field model does not require consideration of unpredictable rock properties or stress inhomogeneities around crack tips. It is shown by benchmarking the models with physical experiments that the numerical assumptions in the phase-field approach do not affect the final model predictions of wing crack nucleation and growth. With this study, we demonstrate that it is feasible to implement the formation of wing cracks in large scale phase-field reservoir models

Auxin export from proximal fruits drives arrest in temporally competent inflorescences

A well-defined set of regulatory pathways control entry into the reproductive phase in flowering plants, but little is known about the mechanistic control of the end-of-flowering despite this being a critical process for optimization of fruit and seed production. Complete fruit removal, or lack of fertile fruit-set, prevents timely inflorescence arrest in Arabidopsis, leading to a previous proposal that a cumulative fruit/seed-derived signal causes simultaneous ‘global proliferative arrest’. Recent studies have suggested that inflorescence arrest involves gene expression changes in the inflorescence meristem that are, at least in part, controlled by the FRUITFULL–APETALA2 pathway; however, there is limited understanding of how this process is coordinated at the whole-plant level. Here, we provide a framework for the communication previously inferred in the global proliferative arrest model. We show that the end-of-flowering in Arabidopsis is not ‘global’ and does not occur synchronously between branches, but rather that the arrest of each inflorescence is a local process, driven by auxin export from fruit proximal to the inflorescence apex. Furthermore, we show that inflorescences are competent for arrest only once they reach a certain developmental age. Understanding the regulation of inflorescence arrest will be of major importance to extending and maximizing crop yields

HDG-NEFEM with Degree Adaptivity for Stokes Flows

This paper presents the first degree adaptive procedure able to directly use the geometry given by a CAD model. The technique uses a hybridisable discontinuous Galerkin discretisation combined with a NURBS-enhanced rationale, completely removing the uncertainty induced by a polynomial approximation of curved boundaries that is common within an isoparametric approach. The technique is compared against two strategies to perform degree adaptivity currently in use. This paper demonstrates, for the first time, that the most extended technique for degree adaptivity can easily lead to a non-reliable error estimator if no communication with CAD software is introduced whereas if the communication with the CAD is done, it results in a substantial computing time. The proposed technique encapsulates the CAD model in the simulation and is able to produce reliable error estimators irrespectively of the initial mesh used to start the adaptive process. Several numerical examples confirm the findings and demonstrate the superiority of the proposed technique. The paper also proposes a novel idea to test the implementation of high-order solvers where different degrees of approximation are used in different elements